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We study the relativistic Hartree approach with the exact treatment of the vacuum polarization in 
the Walecka (a — uj) model. The contribution from the vacuum polarization of nucleon-antinucleon 
field to the source term of the meson fields is evaluated by performing the energy integrals of the 
Dirac Green function along the imaginary axis. With the present method of the vacuum polarization 
in finite system, the total binding energies and charge radii of ie O and 40 Ca can be reproduced. On 
the other hand, the level-splittings in the single-particle level, in particular the spin-orbit splittings, 
are not described nicely because the inclusion of vacuum effect provides a large effective mass with 
small meson fields. We also show that the derivative expansion of the effective action which has 
been used to calculate the vacuum contribution for finite nuclei gives a fairly good approximation. 

PACS number(s): 21 .10.-k,21.60.-n, 13.75.Cs 

I. INTRODUCTION 

It is well known that the relativistic field theory based on the hadrons referred to as quantum hadro- 
dynamics (QHD) has been of great success in describing the ground states of finite nuclei[l|. When the 
energy functional of the relativistic mean field (RMF) is fitted to nuclear saturation, the RMF model 
automatically produces an appropriate order of the spin-orbit splitting of nuclei, the spin-observables of 
the proton-nucleus scattering and energy dependence of the proton-nucleus optical potential. In RMF, 
however, only positive-energy nuclcons arc taken into account in the calculation in spite of the existence 
of solutions with negative energy, which is one of the interesting characters of the relativistic picture. The 
negative-energy states are observables in antinucleonic atom. The bound levels of antinucleonic atom 
are predicted to be much deeper than those of ordinary nucleon since the magnitudes of nucleon-scalar 
and -vector self energies, which are very large, cancel each other to provide the usual binding energy 
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in the nucleon sector, while they add each other in the antinucleon sector (see Ref. |2( and references 
therein). It is important to study if this feature remains to be the case, when we treat the negative-energy 
states explicitly in the QHD framework. In particular, to know how deep the antinucleons are bound in 
the nucleus provides valuable information for the search of the exotic collective antimatter production 
proposed by Grcincr 3] . 

Recently, the gauge invariant nuclear-polarization calculation was carried out in the relativistic 
random-phase approximation (RRPA) based on the RMF theory. It was found that the RRPA eigen- 
statcs with negative energy have a significant role because the transverse form factors to these states 
become considerably large ^ . If antinucleons are deeply bound, the transverse response function has a 
contribution from the antinucleon states with lower energy than twice the nucleon mass. Then, electron- 
scattering or photo-absorption may confirm the large binding of antinucleon. In order to investigate the 
antinucleon state in the QHD model, however, it is required to take into account the effect of vacuum 
polarization of nucleon-antinucleon field to the mean field because of consistency. Namely, the RMF 
model without the vacuum contribution has to be extended to the full one-nucleon-loop approximation, 
which we refer to as the relativistic Hartree approximation (RHA). 

The vacuum contributions and their effects on the bound states of positive energy were investigated 
by several authors within the local-density approximation and the derivative-expansion method 
& y, [l(J, |l]J, [l2( . After refitting the parameters of the model to the properties of spherical nuclei, it was 
found that the RHA calculation can reproduce the experimental data of the binding energies and the 
rms radii of nuclei as well as the RMF can. On the other hand, due to the decrease of the scalar and 
vector potentials by the feedback from the vacuum in the RHA calculation, the effective mass of nucleon 
becomes large and the binding energy of antinucleon becomes small compared with the RMF. 

Despite the finding of the importance of the vacuum effect, however, the exact evaluation of one-loop 
corrections in a finite nuclear system has never been performed. This is an exceedingly difficult task, 
since the exact treatment of vacuum polarization in finite nuclei requires the computation of not only 
the valence nucleon with positive energy but also infinite number of the negative-energy states. In this 



3 3 



the 



context, there is an excellent method developed in the quantum electrodynamics (QED) 
summation of the eigenstates is replaced by the energy integral of the Dirac Green function along the 
imaginary axis in which the Green functions do not oscillate and therefore, it is possible to calculate 
it much faster. Blunden carried out the exact RHA calculation of QHD with 1+1 dimensions and the 
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calculation of quantum solitons model with 3+1 dimensions 



In the present paper we will apply 



this method developed in atomic physics to the RHA calculation of QHD. 

This paper is organized as follows. In Sec. II. A we introduce the effective Lagrangian density used in 
this work and give the outline of the renormalization in the source term of meson fields. The method 
to obtain the renormalized vacuum densities is given in Sec. II. B for baryon density, and in Sec. II. C 
for scalar density. The detail of the computational procedure in the calculation of vacuum correction 
are described in Sec. III. In Sec. IV, the numerical results of the RHA calculation including vacuum 
corrections for baryon and scalar densities are presented for 16 O and 40 Ca and the role of vacuum in the 
properties of nucleus is discussed. In Sec. V, we will also compare our rigorous method with the local- 
density approximation and the derivative expansion. Finally we give summary of the present calculation 
in Sec. VI. 

II. RELATIVISTIC HARTREE APPROACH IN FINITE NUCLEUS 

A. Lagrangian density 

The nucleus is described as a system of Dirac nucleons which interact in a relativistic covariant manner 
through the exchange of several mesons; scalar meson (cr) produces a strong attraction while isoscalar 
vector meson (u) produces a strong repulsion for the nucleon sector. In the present work, we employ the 
Lagrangian density including photon (A) as well as a and oj mesons as 
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-g^Ncr^N ~ g^Nl^^pN - ]^e$ N (l + r^j^A^ipN - SC, (1) 

where 

SC = —CuidftUiv - dyUJ^) 2 + ^(ad^ad^a + a. x a + ^a 2 <r 2 + ^a 3 a 3 + ja 4 a 4 (2) 

denotes counterterms to renormalize the nucleon density which has the divergence from the vacuum. Since 
this Lagrangian density includes the nonlinear coupling terms for a meson, the one-cr-meson loop also 



can contribute to the nucleon density ' 



. In the present paper, however, we neglect this contribution 



for simplification. Assuming that the nuclear ground state is spherically symmetric, the Hartree basis 
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consists of eigenfunctions of the following Dirac equation; 

^jv(r) = 0. (3) 
The meson fields uio(r) and cr(r) satisfy the Klein-Gordon equations, 



J°{tN - 9uU {r) ~ ^e(l + r 3 )A (r)) + ry • V - (mjy + 9aCr(r j) 



(-V 2 + m2) Wo (r) = SArell (r), (4) 
(-V 2 + m>(r) = - 5s /v ron (r) - . 92 ^ 2 (r) - .g 3 ^ 3 (^), (5) 

respectively. The renormalized baryon (p w rcn (r)) and scalar (po- rcn (r)) densities are given by 

/dz 
—. Tr[ 7o G>,r;z)] + (CT) (6) 

Pa re n(r) = J |^Tr[G>, r; z)] + (CT), (7) 

where G H (r, r ; z) is single-particle Green function of the Hartree approximation with the potential terms 
for the Dirac equation Q. The z integrations are carried out along the modified Feynman contour which 

j 2, m| . The divergences arising from these integrals are 



is below the real axis up to nuclear Fermi energy 



removed by the contributions from the counterterms denoted by CT. This procedure will be discussed 
in detail in subsection II. B and II. C. The integral along the Feynman contour can be changed to the 
integral over the imaginary z axis with the additional pole contribution of the positive-energy states up 
to Fermi level. Thus, we may write unrenormalized baryon and scalar densities as 

/ f-Tr[ 7 °G>,r;z)] = ]T ^(r)^(r) - / ^Tr[ 7 °G^(r, r; z)] (8) 
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/ — Tr[G H (r,r;z)] = J2Mr)Mr)+ / — Tr[G^(r, r; z)], (9) 
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respectively. Here, Gy is the vacuum part of the single-particle Green function of the relativistic Hartree 
approximation: 

G H ( ri , r 2 ; z) = Gg(n, r 2 ; z) + G£(n,r 2 ; z) (10) 

F 

Gg(r 1( r 2 ; z) = 2m £ S(z - eiM(ri)&(r a ) (11) 

ei>0 

G*(r 1( r 2 ; z) = V + y ^TiMiM (12) 

The numerical integration for Gy along the imaginary z axis can be carried out straightforwardly because 
there are no poles in the integrand. Although the second terms of right-hand side of © and © have 



divergence, an expansion of the total vacuum correction in the coupling constants g u and g a of meson 
fields verifies that all divergences are contained in the first order of g u for baryon density, and are 
contained up to the third order of g a for scalar density. In the following subsections, we show that these 
divergences are removed by taking the proper counterterms (2) into account. 

B. Vacuum Correction for Baryon Density 

In this subsection, we consider the vacuum correction for baryon density. For the estimation of vacuum 
correction, we will treat proton and neutron on the equal footing to save numerical effort. In order to deal 
with the unrenormalized density containing the divergence, we start with the perturbative expansion by 
uj and <t fields; 

n+ioo t r +ioo 
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where G° denotes the Green function for the free Dirac equation. This expansion allows us to write the 
vacuum polarization as a sum of infinite terms as shown in Fig. 1. 

We note that only the second term in the expansion, whose Feynman diagram is depicted in Fig. 1(c), 
contains an essential divergence of the u> meson self-energy type. This situation is the same as the 
vacuum polarization for electron-positron field in QED except for propagating particle is massive and 
the diagrams including not only vector mesons but also scalar mesons contribute to the correction for 
baryon density. According to Wichmann and Kroll[l^|. the one- loop vacuum correction may be obtained 
by the sum of the finite part of Fig. 1(c) and Fig. 1(a) subtracted by Fig. 1(c). Thus, the renormalized 
baryon density from vacuum, p^f cn (r), is written as 

CW^lW + ^W- (14) 

where p^\ cn (r), which corresponds to the Uehling term[18] in QED, denotes the finite contribution from 
Fig. 1(c) and can be calculated from the renormalized result in the momentum representation. Choosing 
the renormalization point at q 2 = 0, only a wave function conterterm £ w is needed to obtain the finite 



result: 

P&Ur) - / 7^3 ^^oW^lpl 2 f dxx{l - x) In ( 1 + NMl 2 " T) ) , (15) 

where Wo(p) is the Fourier transform of u> meson field CL>o( r )- We estimate the contribution from the 
lowest order of <j£, by this explicit expression numerically. The second term of Eq. I|14|) , which corresponds 



to the Wichmann-Kroll term 
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in QED is the residual finite density expressed by 



P^{r)= P V{r)-p^(r) 

f+ioo j f -+ioo , p 

= J ^- Tr ^° G ^( r ' r > Z)] - j . h Tl[ J ^° G °( r ' x > z)9^°u(x)G°(x, r; z)]. (16) 

In the present work, we evaluate the vacuum correction for baryon density by using this expression 
directly. After the partial-wave expansion by the Dirac angular-momentum quantum number k, each \k\ 
contribution of Eq. (|16|) is still finite. The partial-wave Green function of RHA calculated numerically on 
the imaginary axis is used in the first term, while the analytical form of the partial-wave Green function 
of the free Dirac equation is used in the second term. 

There are several advantages to the method solving the Green function on imaginary axis. For example, 
we can carry out the z integration without taking care of the poles. In addition, we can employ Gaussian 
quadrature for the radial integral in the second term of right-hand side of Eq. (|16f) as well as z integration 
because the Green function on imaginary axis behaves like the modified Bessel function which is not an 
oscillating function. As a result, the vacuum correction can be obtained very fast with high precision 
and it is possible to perform the RHA calculation with practical computational time. 

It should be noted that pffl(r) contain superficially divergent contribution with the three u> mesons 
attached to the baryon loop according to naive power counting, which actually vanishes as a result of 
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current conservation. In QED, how to calculate this term has been discussed by many authors 
and it is well known that this contribution vanishes if the summation over k is restricted to a finite number 
of terms. Then, the k tail contribution is given by the extrapolation. This conclusion is valid for the 
present case since the neutral u> meson couples to the conserved baryon current. Hence Eq. (|14|) together 
with 11511 and (|16f) can be used to calculate the vacuum correction for the baryon density numerically. 

C. Vacuum Correction for Scalar Density 



Next, we consider the vacuum correction for scalar density. The regularization procedure of scalar 
density is performed under the same concept as in the baryon density. One can easily see that the 
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perturbative expansion of the Hartree vacuum correction for scalar density corresponding to Eq. (|13|) 
gives the Feynman diagrams depicted in Fig. 2, where four diagrams from Fig. 2(b) to Fig. 2(e) contain 
divergent contributions. 

The scalar density is renormalized by the counterterms aicr(r) + l/2a 2 cr 2 (r) + l/3a 3 a 3 (r) + 
l/4a 4 o- 4 (r) + l/2C CT 9 M cr(r)^ I CT(r) in SC. The finite part of Fig. 2(c) is calculated fr om the renormalized 
result in the vacuum. The mass (a 2 ) and wave function (Co-) counterterms are required to make the 
contribution finite. The result is given by 



P?ln(r) = I 7^^ r a(p)|% [l|p| 2 - / dx(m% + \p\ 2 x(l - xj) In (l + 

J {2ny [6 Jo V m- N 



(17) 



where cr{p) is the Fourier transform of a(r) and we have chosen q 2 = for the renormalization point. 
The physical vacuum contribution for scalar density is given by the sum of Paren(r) and the residual 
finite density denoted by p a (r): 

P^ e n(r) = P^ln(r) + pi R) (r). (18) 

where the residual finite density p„ -(r) includes the finite contributions arising from Figs. 2(b), 2(d), 
and 2(e). The finite residual vacuum density p&(r) is evaluated by subtracting the contribution of 
Fig. 2(c) and the counterterms from the unrenormalized divergent scalar density: 

/+ioo i 

- / T- Tr t / dxG\r, x; z)g a a{x)G\x, r; z)] 

J-ioc 2tti J 

-(ai + a 3 a 2 {r) + a 4 a 3 (r))/ gr7 , (19) 

where the coefficients ai, «3, and of the counterterms are given by 

'•+ 100 dz 
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27r .Tr[G°(r,r;z)} (20) 



/ dz f 
—Tr[j dx 1 dx 2 G (r,x 1 ;z)G (x 1 ,x 2] z)G (x 2 ,r;z)} (21) 

f +io ° dz f 

/ — Tr[ / dx 1 dx 2 dx 3 G°(r,x 1 ;z)G (x 1 ,x 2 ;z)G (x 2 ,x 3 ;z)G°(x s ,r;z)}. (22) 
J_ ioo 2tti J 



We note that a±, a 3 , and are independent of r though it remains in the right-hand side of these ex- 
pressions. Performing the partial-wave expansion of free Green functions, however, each \k\ contribution 
to the coefficients depends on the radial part r. With these partial-wave subtraction terms, the scalar 
density is calculated for each The net effect of scalar density generated from the vacuum polarization 
are obtained by taking the extrapolation K maa: | — > oo. 
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III. COMPUTATIONAL DETAIL FOR VACUUM CORRECTION 

As has been explained in the previous section, pi^rcn and p^ 1 ren contributions are calculated from the 
explicit forms (|15|l and l|17f) . They contribute by a small amount as discussed below. Therefore, here we 
explain mainly the estimation of the residual contributions defined by Eqs. (|l(jfl and (|19|) . We write the 
radial part of the residual baryon and scalar densities as 

W = "£%-en,|«|« (23) 
1*1=1 



P R 

ruj rcn 



and 

^renW= £ ^rc„, W M> (24) 
1*1 = 1 

respectively. Here, pff ren | K | (p^ron |«|) represents the contribution from |k| = ±k to net residual baryon 
(scalar) density. We compute p^ Tcn i re i and p^\. cn using the angular momentum decomposed expression 
of Eqs. i|ltjfl and i|19[l . The detail of our calculation is the following. For the respective contributions 
of \k\ in Eqs. Ijlbj) and l|19|l . we carry out the integral over iz using Gaussian quadrature. The radial 
Green functions with the imaginary energy are given in terms of the two solutions of the Dirac equation; 
the solution regular at r = and the solution regular at r — oo. Several values of the upper and lower 
limits of the integral over z are selected depending on the radius parameter r around 20i ~ 50i GeV. 
The vacuum densities p^ ron | K |( r ) ano - Paicn K|( r ) ^ or z max oo are extrapolated from the resulting 
integrated values. Then, the sum over k| is performed until the cutoff |K maa; |. The sum over |k| does not 
converge so rapidly. In the present work, |K ma :r| in baryon density is 33 for both of 16 and 40 Ca, while 
l^maa;! hi scalar density is 36 for 16 and 48 for 40 Ca. Finally, we have extrapolated \k\ contribution to 
larger values to obtain the convergent vacuum densities. As an example, the contributions from |k| =5, 
10, 15, 20, 25, and 30 for in 16 are shown in Fig. 3. 

An accuracy test for the computed vacuum correction of the baryon density p^ rea (r) is provided by 
the requirement that for the total induced vacuum correction should vanish due to the conservation of 
baryon density, 

/>OC 

AB = 4n drr 2 pK ICU (r) = 0. (25) 



In the present calculation with various QHD parameters, we found AB ^10 3 as a typical values both 
for 16 O and 40 Ca. For the scalar density, on the other hand, we have no requirement of conservation 
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laws. However, it would be reasonable to expect the numerical error in scalar density is of the same 
order of magnitude as the one in the baryon density. 

In Sec. II. B, we mentioned that the divergence of unrenormalized baryon density in the present model 
has the same structure as unrenormalized vacuum charge density of electron-positron field in the QED 
correction. However, the dependence on the partial-wave contribution is largely different from the QED 
case: a large |K maa; | is required to achieve the convergence as seen in Fig. 3, while for the renormalized 
charge density in QED only the term with |K ma a;| = 1 is a good approximation 13J. In addition, it should 
be pointed out that Eqs. I|15l) and l|17|) are negligible in the present calculation, namely, pu\cn and piHcn 
are an order of magnitude smaller than the residual densities in contrast with the QED calculation, 
in which the corresponding contribution known as the Uehling effect has a dominant contribution^]]. 
These differences from QED seem to be caused by the u> — a couplings as well as the large coupling 
constant of the nuclear force. In particular, the baryon density, the source of the oj meson, is strongly 
influenced by the a meson. This is shown in Fig. 4 where the baryon densities p£ renM with and without 
a meson are compared. We can see there that the vacuum correction for baryon density is negligibly 
small if the lo meson only is taken into account. We found that the baryon density induced by the 
vacuum polarization comes from the diagrams with a meson self-energy insertions mainly. 



Here we show the results of the relativistic Hartree calculation with a rigorous treatment of renormal- 
ized densities in 16 O and 40 Ca based on the Lagrangian density of Eq. Q). The numerical procedure of 
the present RHA calculation is similar to that used in the conventional RMF calculation: first, the Dirac 
equation @ for only the valence nucleons is solved under the external w and a fields. Second, using the 
same potential, we calculate vacuum densities by means of the Green function described in the previous 
section. Third, the equations of motion of mesons, © and Q, are solved with the source terms due 
to the valence nucleons and vacuum contributions. Substituting these results in the Dirac equation, we 
complete one iteration step. Using the same method as indicated in the previous section, the energy 
shift due to the vacuum polarization is estimated by 



IV. RESULTS AND DISCUSSION 




(26) 



CT = -l/4C(Vw ) 2 + l/2C CT (Va) 2 + a x a + l/2a 2 cj 2 + l/3a 3 a 3 + l/4a 4 CT 4 , 



(27) 
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is calculated at each step of iteration. The iteration is continued until the total binding energy of the 
nucleus -Etotai = E mcscm + E va \ cncc + Eyp converges showing self consistency. 

The QHD parameter set is chosen so as to reproduce reasonably well the experimental values of the 
total binding energies, the rms radii, and the single-particle energies for both of 16 and 40 Ca. In 
the second column of Table I, we give the results with the coupling constants and masses g a = 7.38, 
52 = 7.90, g 3 = 3.20, and m a = 458.0 MeV for a meson, and g u = 9.18 and m u = 783.0 MeV for u> 
meson. The results of the RMF calculation with the parameter set TM2 2(| and experimental data taken 
from Ref. | are also shown in the third and last columns, respectively. We see that our total binding 
energies including the vacuum correction and rms radii are similar to those of RMF and agree with the 
experimental data well. However, it should be noted that the present result is produced by using smaller 
coupling constants than those of RMF. Hence, the scalar and vector potentials in the present RHA differ 
considerably from those of RMF. The results of the present RHA in 16 and 40 Ca are compared with 
the results of RMF in Figs. 5 and 6 for scalar and vector potentials, respectively. One can see that the 
holds in the RHA are considerably small compared with those in the RMF. 

The vacuum contribution plays a crucial role for the generation of the weak meson fields. In order to see 
this, we plot cr-meson field in nuclear matter as a function of the coupling constant g a while keeping the 
other parameters fixed in Fig. 7. The figure shows that the vacuum correction contributes destructively 
to the valence contribution and it is impossible to obtain a strong meson field unless a very large coupling 
constant is used. Of course, we could choose a parameter set with large coupling constants. However, 
such a parameter set would have produced an unstable solution in the RHA calculation, because if the 
deeply-bound antinucleon states are produced by the strong to and a fields with large coupling constants, 
then the vacuum effect, which works in the opposite direction, also becomes large to suppress the strong 
fields. Such a RHA solution is not realistic because it is unstable even for the trivial fluctuation of the 
density. Thus, we have to choose a parameter set which produces a weak field in the self consistent 
iteration. 

The RMF reproduces reasonably well the observed tendency of the single-particle spectra due to the 
small effective mass of the nucleon, m* N {r) — mjf + g a a{r). The fact that the scalar field is suppressed 
in RHA results in the large effective mass. As a result, it raises a problem in fitting the single-particle 
energies. As seen in Table I, the energy splittings in the single-particle states of the present RHA are very 
small and they are unlikely to agree with the experimental values. This problem has been already known 
by RHA calculation using the local-density approximation and the derivative expansion to estimate the 
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vacuum correction|8j,lllJ. It is difficult to resolve this problem by RHA with ordinary QHD models used 
up to now. 

Thus, the QHD models require a mechanism for the spin-orbit splittings other than the small effective 



mass. One suggestion is made in Refs. 



12, 



where a tensor-coupling of u meson is introduced in order 
to provide the spin-orbit splittings. However, the omega meson coupling with the nucleon is known to 
be dominated by the vector coupling in the nucleon-nucleon potential. Another candidate to solve this 
problem may be the possibility of the finite pion mean field in the relativistic Hartree framework, which 
was suggested to provide the spin up and spin down partners to have large energy separations 22J. It is 
interesting to extend the present RHA calculation with these effects taken into account. This is certainly 
a subject to be worked out in future study. 

V. COMPARISON WITH THE PREVIOUS METHOD FOR VACUUM POLARIZATION 



f the negative-energy nucleons for finite nuclei was first estimated by the local-density 

ion .... nnff 

■i HI- l| • It was developed further by applying the derivative-expansion method [y, |g, 1111 



The effect of thejregative-energy nucleons for finite nuclei was first estimated by the local-densit) 
approximation 

12j ■ In this section, we compare our results densities induced by the vacuum polarization with those of 
the local-density approximation and the derivative-expansion method. The local-density approximation 
uses as input the result of the infinite nuclear matter. In this approximation, the vacuum correction is 
given by 

P*re ( n DA) (r) = ~ [m# \n(m N /m N ) + l/3m% - 3/2m 2 N m N + 3m N m% - ll/6m*J] , (28) 

7T 

and decreases the scalar density in nuclear interior. The vacuum does not change the baryon density 
because of the conservation of the baryon number. In the derivative-expansion method, on the other 
hand, the presence of the derivative term allows non-vanishing correction for baryon density as well as 
for scalar density: 

pZ P r£ E \r) = * V • In (^1) Vu (r), (29) 
P^ E) (r) = P v a P r ^ DA \r) - • In Va(r) 

Z7T \ WIN J 

- A 2 9 \ ( , (^(0) 2 + fi 2 9 \, , (Vc(r)) 2 , (30) 
4ir z m N (r) bir z m N (r) 

where the leading order of the derivative terms is taken into account. The baryon and scalar densities 
induced by the vacuum polarization are given in Fig. 8, together with those from the local-density 
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approximation and the derivative expansion. For purpose of comparison, we assume there the same 
potential in evaluating the vacuum polarization. We can see that the densities obtained by the local- 
density approximation are corrected significantly not only for baryon density, which vanishes in this 
approximation, but also for scalar density. The both of scalar and baryon density-profiles obtained by 
the present calculation are in a surprisingly good agreement with those of the derivative expansion. 

However, this cannot be always the case|l(j and it is possible to attribute the excellent agreement 
between our method and the leading-order derivative expansion to the character of a— u) model. Consider, 
for example, the vacuum correction in the baryon density without a meson. In this situation, the vacuum 
correction using the present method can be significant with large coupling constant of u> meson. As found 
in Eq. (|29|l . on the other hand, the vacuum correction from the derivative expansion vanishes exactly for 
m* N (r) -> mjv- Hence, we find that a meson plays important role in the agreement between our method 
and the leading-order derivative expansion. Thus, the present calculation supports that leading-order 
derivative expansion is greatly useful for the estimation of the vacuum correction in RHA. 

In the present calculation of the full RHA for finite nuclei, the vacuum-polarization corrections (|15|) 
and (|17|l to the meson propagators are implicitly taken into account to all orders through iteration to 
achieve the self-consistency in the relativistic Hartree approximation. There the unphysical pole in the 
meson propagator at finite momentum transfer, known as Landau ghost, may affect the present numerical 
resultsj^, l24| through integration over the momentum transfer. However, this unphysical effect is not 
significant in the RHA of finite nuclei, because even if the H15[) and l|17f) terms are totally ignored in the 
calculation, the final results of scalar and vector densities do not change appreciably. The good agreement 
with the results of the derivative expansion, where the Landau ghost plays no role, also implies that the 
unphysical effect is negligible. 

VI. SUMMARY 

We have developed a rigorous method to calculate vacuum-polarization effects in relativistic Hartree 
approach. The renormalized baryon and scalar densities have been evaluated within a practical compu- 
tational time replacing the summation of the Hartree basis by the numerical integral of the Dirac Green 
function over the imaginary energy. We have obtained numerical results that the vacuum corrections 
for baryon and scalar densities are non-negligible in the RHA calculation. In particular, we have seen 
that the vector and scalar potentials in the RHA are largely suppressed by the feedback of the vacuum 
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polarization to the mean fields. Such suppressed potentials of the u> and a fields affect the theoretical 
calculations for binding energies of antinucleon and the various sum rules. 

Our results with the Walecka model have reproduced the experimental binding energies and rms radii 
of 16 and 40 Ca nicely after adjusting the parameters. However, it was impossible to find a QHD 
parameter set to reproduce spin-orbit splittings in accordance with the observed data and required by 
nuclear shell model. In the a — u> model, the main attraction is caused by the large a mean field, which 
provides small nucleon effective mass in the finite nuclei. However, the negative-energy nucleons do not 
like to change its mass from the free value. Hence, the negative-energy nucleons try to keep the nucleon 
mass at the free value. For the net effect, the effective nucleon mass remains to be quite large. Once the 
effective nucleon mass becomes large, the spin-orbit splittings in the single particle spectra come out to 
be very small. The QHD type effective theory based on the a — lo mesons need to include new type of 
interaction terms and/or to go beyond the RHA approximation to solve this problem. 

We have found that our result of RHA calculation is very similar to that in Refs. (], 11 1 where the 
derivative-expansion method was used to estimate the vacuum polarization. In particular, it has been 
shown that the agreement of density-profiles of the vacuum correction is quite good. Thus, the validity 
of this approximation has been confirmed by the present calculation. 
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TABLE I. The total binding energies, the rms charge radii, and the single-particle energies in 16 and 
40 Ca. 
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FIG. 1: Graphical representation of the baryon density in the self-consistent relativistic Hartree approximation. 
The double and single lines denote the Hartree propagator and the free propagator of nucleon, respectively. The 
wavy and dotted lines with cross represnt the vector and scalar potentials given by the previous step of Hartree 
iteration respectively. The divergence contained in (a) is caused by the contribution from (c). 
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FIG. 2: Graphical representation of the scalar density in the self-consistent relativistic Hartree approximation. 
Same notations as in Fig. 1 are used. The divergence contained in (a) is caused by the contributions from (b) to 
(e). 





FIG. 3: The contributions from |k| =5, 10, 15, 20, 25, and 30 to (a) the baryon density rm ,\ K \ (r) and (b) the 
scalar density p? relliW (r). 
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FIG. 4: Comparison between the baryon densities with and without a meson contribution. The latter is given 
by solving the Dirac Green function with the u meson field only. 
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FIG. 5: The scalar potentials for 16 O in the left panel and for 40 Ca in the right panel. The scalar potential with 
the vacuum polarization is shown by solid curve, while the one without it by dashed curve. 



21 




FIG. 6: The vector potentials for O in the left panel and for Ca in the right panel. The scalar potential with 
the vacuum polarization is shown by solid curve, while the one without it by dashed curve. 
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FIG. 7: The scalar potential in nuclear matter, cr meson mass m CT = 458.0 and Fermi momentum &f = 1-42 are 
employed, er denotes the ordinary a meson field generated from the valence nucleons while avac denotes the 
contribution from the vacuum. Due to the cancellation between them, the net a meson field does not increase 
smoothly with the coupling constant, g a . 



